##Exerccios da Aula 7


#1) H uma crena do senso-comum de que basta duplicar a altura aos dois anos de idade para descobrir a altura que uma criana ter 
#quando se tornar adulta. Imagine uma amostra de oito pessoas adultas, tomada ao acaso de uma populao, para as quais temos as 

#Altura aos dois anos (in) 39 30 32 34 35 36 36 30 
#Altura adulto (in) 71 63 63 67 68 68 70 64 


bebe=c(39, 30, 32, 34, 35, 36, 36, 30) 
adulto=c(71, 63, 63, 67, 68, 68, 70, 64)

#1.Faa uma regresso da altura na idade adulta em funo da altura na idade de dois anos. H uma relao significativa? Justifique.


###lm- ajusta modelo linear gausiano##
alt= lm(adulto~ bebe)
anova(alt)#0.00386 #modelo alt x modelo nulo

#2.Faa um grfico de disperso com os dados, e acrescente a linha da regresso, e a linha esperada pela crena.

plot(adulto~bebe)
abline(alt, col="red")

crena= bebe*2
crena


alt2=lm(crena~bebe)# modelo alt X modelo da crena
abline(alt2, col="blue") 
# de acordo com a crena os bebes seriam bem maiores !

#3.Calcule os intervalos de confiana dos coeficientes.
help(confint)# computa o intervalo de confiana de um ou mais modelos
help(coef)# extrai m coeficiente do modelo

confint(alt)
confint(alt2)

#4.Seus resultados corroboram a hiptese do senso comum? Por que?
# no, o esperado pela crena  diferente do observado 


###Seriemas e Carcars###


#Use o conjunto de dados Aves no Cerrado para avaliar se o nmero de avistamentos de seriemas  afetado pelo nmero de Carcars avistados, em cada 
#fisionomia de cerrado. Para cada pergunta abaixo indique os comandos que usou para respond-la (alm de respond-la  ). 

av=file.choose()
aves=read.table(av,header=TRUE, sep=";",as.is=TRUE)
head(aves)#ok

#1.H relao entre o nmero de avistamentos das aves em alguma das fisionomias amostradas?
# Fisionomias = Ce, CC, CL/ Aves= urubu, carcara, seriema

#CE- urubu~carcara, urubu~seriema, carcara~seriema


ce1=lm(aves$urubu~aves$carcara,subset=aves$fisionomia=="Ce")
anova(ce1)#0.391
ce2=lm(aves$urubu~aves$seriema,subset=aves$fisionomia=="Ce")
anova(ce2)#0.432
ce3=lm(aves$carcara~aves$seriema,subset=aves$fisionomia=="Ce")
anova(ce3)#0.4803

#CC
cc1=lm(aves$urubu~aves$carcara,subset=aves$fisionomia=="CC")
anova(cc1)#0.977
cc2=lm(aves$urubu~aves$seriema,subset=aves$fisionomia=="CC")
anova(cc2)#0.1640
cc3=lm(aves$carcara~aves$seriema,subset=aves$fisionomia=="CC")
anova(cc3)#0.5426

##CL
cl1=lm(aves$urubu~aves$carcara,subset=aves$fisionomia=="CL")
anova(cl1)#0.5788
cl2=lm(aves$urubu~aves$seriema,subset=aves$fisionomia=="CL")
anova(cl2)#0.002305 ############ Significativo !!!
cl3=lm(aves$carcara~aves$seriema,subset=aves$fisionomia=="CL")
anova(cl3)#0.01853

##s h relao significativa entre  urubu e seriema na fisionomia CL 

#2.H diferenas na relao entre o nmero de avistamentos de seriemas e carcars entre as fisionomias do Cerrado?

plot(aves$seriema[aves$fisionomia=="Ce"]~aves$carcara[aves$fisionomia=="Ce"])
pp1=(aves$seriema[aves$fisionomia=="Ce"]~aves$carcara[aves$fisionomia=="Ce"])

p1=lm(aves$seriema[aves$fisionomia=="Ce"]~aves$carcara[aves$fisionomia=="Ce"])
abline(p1)
p2=lm(aves$seriema[aves$fisionomia=="CC"]~aves$carcara[aves$fisionomia=="CC"])
pp2=lm(aves$seriema[aves$fisionomia=="CC"]~aves$carcara[aves$fisionomia=="CC"]
abline(p2, col="red")
p3=lm(aves$seriema[aves$fisionomia=="CL"]~aves$carcara[aves$fisionomia=="CL"])
abline(p3, col="blue")

# aparentemente CL tem uma relao diferente de CC e Ce, o nmero de espcies 
#carcara cai muito mais drasticamente com o aumento do nmero de seriemas avistadas


##Resduos de Iris##


#Uma maneira simples de descontar o efeito indesejado de uma co-varivel  usar os resduos de uma regresso de sua varivel de interesse em funo desta co-varivel.
# Use as observaes da espcie Iris setosa, no objeto de dados iris para: 

iris
head(iris)

#1.Ajustar e interpretar uma regresso linear da largura da spala em funo do comprimento da spala.

# comprimento = Length, largura=width

plot(Sepal.Length~Sepal.Width , data=iris,subset=Species=="setosa")
sep= lm(Sepal.Length~Sepal.Width, data=iris,subset=Species=="setosa")
abline(sep, col="red")
# aparentemente se ajusta bem
anova(sep)#p=6.71e-10 ***-  significativo, a largura da slapa  melhor explicada pelo comprimento da sep do que pelo modelo nulo de que no h relao
summary (sep)# r^2 = 0.5514- mostrando que 55% da largura da spala  explicada pelo comprimento da mesma


#2.Ajustar e interpretar a mesma regresso, mas descontando do efeito do comprimento da ptala de cada varivel.

desc.p1=lm(Sepal.Width~Petal.Length,data=iris,subset=Species=="setosa")
desc.p2=lm(Sepal.Length ~Petal.Length,data=iris,subset=Species=="setosa")

##resduos - 2 local de names
names(desc.p2)
resi2=desc.p1[2]
resi2# res do comp

head(iris)
summary(iris)

#transformando em numerrico
res.comp=residuals(resi2)
class(res.comp)
Sepal.Length


plot(Sepal.Length~res.comp, data=iris,subset=Species=="setosa")
efeito1= lm(Sepal.res.comp~resi2, data=iris,subset=Species=="setosa")
plot(Sepal.Width~resi2, data=iris,subset=Species=="setosa")
efeito= lm(Sepal.Length~resi2, data=iris,subset=Species=="setosa")




